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ABSTRACT 



We study the distribution of orbital eccentricities of stars in thick disks generated by the heating of a pre-existing thin stellar disk 
through a minor merger (mass ratio 1:10), using N-body/SPH numerical simulations of interactions that span a range of gas fractions 
in the primary disk and initial orbital configurations. The resulting eccentricity distributions have an approximately triangular shape, 
with a peak at 0.2-0.35, and a relatively smooth decline towards higher values. Stars originally in the satellite galaxy tend to have 
higher eccentricities (on average from e = 0.45 to e = 0.75), which is in general agreement with the models of Sales and collaborators, 
although in detail we find fewer stars with extreme values and no evidence of their secondary peak around e = 0.8. The absence of 
this high-eccentricity feature results in a distribution that qualitatively matches the observations. Moreover, the increase in the orbital 
eccentricities of stars in the solar neighborhood with vertical distance from the Galactic mid-plane recently found by Diericxk and 
collaborators can be qualitatively reproduced by our models, but only if the satellite is accreted onto a direct orbit. We thus speculate 
that if minor mergers were the dominant means of formating the Milky Way thick disk, the primary mechanism should be merging 
with satellite(s) on direct orbits. 
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1. Introduction 



ISales et all (l2009t) proposed using the orbital eccentricity of stars 
in a thick disk as a way of constraining the main physical mech- 
anism responsible for its formation. They pointed out the dis- 
tinctive differences between the eccentricity distributions pre- 
dicted for a comprehensive range of formation mechanisms for 
the thick disk; the four mechanisms they considered were: the 
heating of a pre-existing thin disk by a varying gravitational 
potential in the thin disk, whi ch causes radial mig ration and a 

l2008ah: heating of 



thickening of the stellar disk (Roskaret al 



Villalobos & Helmil 
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a pre -existing thin disk by minor mergers 
2008); accretion of disrupted s atellites dAbadi et al.H2003l) ; and 
formation by gas-rich mergers dBrook et al.ll2004| 120051) . Their 
results can be summarized as: (1) the accretion scenario predicts 
a broad distribution of stellar eccentricities, e, which are quite 
symmetric with respect to its peak at around e = 0.5; (2) the ra- 
dial migration scenario produces a thick disk whose stars have 
eccentricities peaking around e = 0.25 - 0.3, and whose distri- 
bution is quite narrow and symmetric around its peak; (3) the 
heating of a pre-existing thin disk by a 1:5 mass ratio merger 
produces a distribution with a peak around e = 0.25, with a sec- 
ondary peak and tail at high eccentricities composed mostly of 
stars originally in the satellite galaxy; and (4) gas-rich mergers 
produce a result quite similar to the previous scenario, except for 
the absence of the secondary peak. The strength of this work is 
that it considers a number of scenarios for the formation of the 
thick disk - each of which apparently produces a unique (and 
testable) signature. 



A comparison of the eccentricity distributions predicted by 
these four models with that of stars populating the thick disk 
in the solar neighborhood led iDierickx et al.l d2010l) to reject 
the above-mentioned models (1), (2) and (3), which leads them 
to propose gas-rich mergers at early epochs as the most prob- 
able mec hanism for th e form ation of the Galactic thick disk. 
Similarly, Wil son et al.l d2010h found that the eccentricity dis- 
tribution of a sample of solar neighborhood thick disk stars from 
RAVE survey data appears to be most consistent with the gas- 
rich merger scenario. Within the limitations of any of these com- 
parisons, their conclusions appear well justified. However, the 
question arises of whether this is the complete story. 

To make this comparison based on a comprehensive analysis 
of the effects of minor mergers, a wide range of initial orbital 
parameters needs to be considered in numerical simulations. The 
aim of this Letter is to compare the observed distribution of the 
eccentricities of thick disk stars in the solar neighborhood with 
those predicted by twelve 1:10 minor mergers simulations, with 
and without gas in the primary, i.e. the most massive galaxy, and 
in the satellite, that have a range of initial orbital characteristics. 

2. Models 

The simulations we analyzed are part of a set of minor merger 
simulations with a mass ratio of 1:10, described extensively in 
IChilingarian et alj d2010l) . The primary and the satellite galaxy 
each consist of a spherical non-rotating halo and a central bulge, 
both modelled as Plummer spheres, and a stellar and an optional 
gas disk represented by Myamoto-Nagai density profiles. Both 
the primary and the satellite galaxies span a range of morpholo- 
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gies, with bulge-to-disk mass ratios of between 0.2 and 0.25 
and gas-to-stellar disk fractions, f gas , from to 0.2. We ana- 
lyzed twel ve of the encounte r s label led gSOdSO, gSadSa, and 
gSbdSb in IChilingarian et al.1 d2010h containing, respectively, 
f g as=0, 0.1, and 0.2. A total of JVtot = 528 000 particles were 
used in all simulations, distributed between the primary galaxy 
and the satellite. We tested the dependence of the results on the 
number of particles used by running an additional simulation 
with 2.5 times more particles, A^tot = 1 320 000. After the 
particles in each simulation had relaxed, the two galaxies were 
placed at an initial distance of 100 kpc, with a variety of ini- 
tia l relative velocitie s (id=0 1dir, Olret, 02dir, 02ret in Table 9 
of IChilingarian et all (120101) ). to simulate different orbits. For a 
complete description of the initial galaxy models, the number of 
particles employed, and the orb ital characteristics, we re fer the 
reader to Tables 1, 3, 6, and 9 in lChilingarian et all d2010l) . 

All simulations were run using the Tree-SPH code described 
in lSemelin & Combesl (12002). The gas is assumed to be isother- 
mal at a temperature of 10 4 K. Prescriptions for star formation 
and feedback from supernovae explosions are also included. The 
equations of motion are integrated using a leap-frog algorithm, 
with a fixed time step of 0.5 Myr. For the evaluation of the gravi- 
tational forces, a softening length e=200 pc (or e=150 pc for the 
high-resolution simulation) is employed. With these choices, the 
relative error in the conservation of the total energy is about 10~ 6 
per time step. 

3. Results and discussion 

We analyzed the resulting eccentricities from the orbits of the 
stellar particles in our modeled galaxies during the period 1-2 
Gyr after the completion of the merger. This is sufficient for the 
orbits to be relatively stable. For each particle, we computed the 
minimal (r m ,„) and maximal distance (r max ) from the galaxy cen- 
ter during this period of time, and defi ned the eccentricit y as 
e = (r max - r min )l(r max + r min ). As done in lSales et aD (120091) . we 
restricted our analysis to a cylindrical region around the galaxy 
center, at 2 < R/R d < 3 and 1 <| z/zo l< 3, where R/R d 
and z/zo are the cylindrical coordinates of a star relative to the 
disk scale-length and the scale height of the thick disk. The disk 
scale-length varies slightly from one model to another, depend- 
ing mostly on the gas fraction in the primary disk, ranging from 
Rd = 4.3 to R c / = 5.4 kpc. The average value of the thick disk 
scale-height in the volume under consideration is zo = 2 kpc 
(fitted with an exponential), whereas the exact value depends on 
both the gas fraction of t he primary and the initial orbital config- 
uration (lOu et al.ll20"T0ah . 

For a representative set of models, the resulting distributions 
of orbital eccentricities of the old stars^in the remnant galaxy are 
shown in FigJT] Independent of the exact properties of the galax- 
ies in the simulations and of the initial orbital configurations, a 
number of common features emerge from the simulations: 

- the overall distribution of eccentricities shows a peak at low 
values, between e = 0.2 - 0.35, followed by a relatively 
smooth decrease towards higher values; 

- stars initially in the satellite galaxy span a wide range of ec- 
centricities, from e = 0.45 to e — 0.75, with the exact values 
depending on the orbital characteristics and gas fractions; 

1 By "old stars", we mean stars already present in the galaxy before 
the interaction, whereas by "new" stars, we mean stars that are formed 
during the interaction from the gas, if initially present in the disk. New 
stars are distributed in both a thin disk (scale height of hundreds of 
parsecs) and the bulge. They do not contribute to the thick disk. 



- we never find a secondary peak at high e ccentricities, simil ar 
to that between 0.7 and 0.95 reported in lSales et aD (120091) . 

Our results ar e quite different from the distribution found by 
ISales et all (120091) for minor mergers (labelled "heating" in their 
Figure 3) - they resemble more what they found for a thick disk 
formed in situ by gas-rich mergers early in the evolution of the 
disk (labelled "merger" in their paper). What is the cause of this 
difference from our simulations? 

In principle, many parameters can contribute to this differ- 
ence. The mass ratio of the encou nters is different, f or example: 
we analyze 1:10 mergers, whereas lSales et aD (120091) considered 
a 1:5 ratio. A relatively more massive satellite would lead to a 
greater fractional contribution of its stars to the total mass in 
the thick disk at a given height above the galaxy mid-plane. To 
investigate the impact of a more massive accretion event, we an- 
alyzed the case of a primary galaxy accreting consecutively two 
satellites^, each having a mass that is one tenth of the mass of 
the primary galaxy. In this case, the total mass acquired through 
the two minor mergers is equivalent to that acquired in a single 
merger of mass 1:5. Using the two consecutive mergers with a 
mass ratio of 1:10 each to investigate the impact of a 1:5 minor 
merger is appropriate because it has bee n shown that the effect 
of multiple min or mergers is cumulative dBournaud et al.|[2007b 
lOu et al.ll2010al) . This being the case, we analyzed several pos- 
sible consecutive merger configurations. We considered mergers 
of two dSO galaxies on a gSO, mergers of two dSa galaxies on a 
gSa, and mergers of two dSb galaxies on a gSb; simulating both 
direct and retrograde encounters. The result of this analysis is 
summarized in Fig. 2, which represents the eccentricity distribu- 
tion of stars that results from a 2x(l:10) merger of a gSb galaxy 
and two dSbs on direct orbits. We find that the consecutive minor 
mergers cause an increase in the contribution of satellite stars to 
the high eccentricity tail of the distribution. We note also that the 
overall distribution tends to shift towards higher values of the ec- 
centricity (this is clear when comparing the results of the direct 
1:10 merger for f gas = 20% in Fig. 1 with Fig. 2). In the six cases 
we analyzed, we never found the secondary peak identified by 
Sales et al. 2009. The mass of the satellite certainly tends to in- 
crease the contribution of satellite stars to the high eccentricity 
tail of the distribution, but it appears insufficient in itself to cre- 
ate a second peak. Other parameters must play a role, perhaps a 
decisive role in creating a second peak. For example, the orbital 
configuration of course has an effect, but this difference may also 
depend on the details of the minor merger model analyzed by 
Sales et al. (2009) which consists of two galaxies (a primary and 
its satellite) embedded in dark halos with a NFW profile. A dif- 
ferent central concentration in the dark matter profile may cause 
dynamical friction to have a slower/delayed impact allowing the 
satellite to survive longer, thus the orbit to become more eccen- 
tric for a more massive and/or centrally concentrated satellite. 
All of this may result in a final eccentricity distribution of stars 
in the remnant galaxy that is different from and extends to higher 
values than those we found. 

We do not mean to claim that the distributions of stellar ec- 
centricities shown in Fig[T]represent the general outcome for mi- 
nor mergers, only that they are a possible outcome, compatible 
with the eccentricity of thick disk stars in the solar neighborhood 
(iDierickx et al.ll2010t IWilson et al"]l2010t) . Nevertheless, the ec- 
centricity distribution found in our simulations, for the direct 
merger with a gas fraction of 20% in the disk, is r emarkably sim- 
ilar to those observed in the Milky Way thick disk. lDierickx et all 

2 Hereafter, we denote simulations of these consecutive mergers as 
"2x(l:10)". 
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Fig. 1. Eccentricity distribution of stellar particles in a number of 1 : 10 minor mergers within a cylindrical shell at disk scale lengths, 
Rd, from 2 < R/R c / < 3, and at heights, zo, from 1 < z/zo ^ 3. From left to right: The gas-to-stellar disk mass fraction increases 
from f g as=0, 0.1, and 0.2 in each panel. Mergers on direct orbits are shown in the top row, ones on retrograde orbits are shown 
in the panels of the bottom row. We indicated both the contribution of stars initially in the satellite galaxy (red histogram) and 
the distribution of all the stars (grey histograms). The dashed histogram in the upper-left panel corresponds to a high-resolution 
simulation, employing a total of 1 320 000 particles and a softening length e=150 pc. The results of the higher resolution simulation 
are in qualitatively in agreement with those of the lower resolutions simulations. 



( 120101) and lWilsonetal] d2010l) found a peak at low eccentrici- 
ties with a smooth decline at higher values as we have in our 
simulations (Fig. [3]). 

To investigate the properties of thick disks formed by minor 
mergers, we have measured the distribution of eccentricities of 
stars at different heights from the disk mid-plane. This analysis 
was motivated by observations of changes in the eccentricity dis- 
tributions of stars in the solar ne ighborhood with distance from 
the plane dDierickx et al.ll20loh . We restricted our analysis to 
stars in the region between 2 < R/Rj < 3 (FiglUi. We have found 
that the value of the peak in the overall eccentricity distribution 
decreases with increasing gas fraction of the primary disk (cf. 
FigfJJand FiglUi. The gravitational force exerted by a gas com- 
ponent can indeed str ongly affect the verti cal stellar distrib ution 
dBanerjee et~al1l2008l) . limit disk heating (lOu et alJl2010ah. and 
affect the kinematics of the stars in the disk dOu et alj|2010bl) . 



Interestingly, we also find that the eccentricity distributions 
as a function of z show a dependency on the initial orbit of the 
satellite galaxy. If the satellite is accreted along a direct orbit (i.e. 
if the z component of the orbital angular momentum is parallel to 
the primary internal angular momentum), the eccentricities show 
a trend with the vertical distance from the mid-plane: the larger 
the distance from the mid-plane, the greater the typical stellar 
eccentricity. This trend can be explained in terms of a rotational 
lag of the thick di sk as a function of z, whi ch can be induced 
by minor mergers (IVillalobos & Helrmll2008l) but who se ampli- 
tude an d behavior strongly depend on the satellite orbit dOu et alj 
1201 Obi) . In direct encounters, the strong coupling of the satellite 



motion to the disk internal rotation induces a higher perturba- 
tion of the kinematics of stars in the primary disk, increasing 
the overall radial velocity dispersion and decreasing the overall 
rotational support. This effect also manifests itself in the ampli- 
tude of the rotational lag of the thick disk. As the stars in the 
thick disk mostly originated in the thin disk of the primary, the 
lag must increase with z- For retrograde orbits, however, the lack 
of strong gravitational coupling between satellite and disk pro- 
duces a much smaller effect - stars in the thick disk rotate slower 
than those in the thin disk, but this lag does not show any depen- 
dence on vertical height, except at very large mid-plane distances 
(above z > 5zo). The absence of a trend between the rotational 
support and z for retrograde encounters is reflected directly in 
the distribution of eccentricities, which does not change signifi- 
cantly with increasing z- In particular, the peak of the distribution 
is constant at | z/zo |<5, and increases only at even larger disk 

heights. 

Thus the observational trend found bv lDierickx et al.l d2010l) . 
where stellar orbits become more eccentric as one moves away 
from the Galactic plane, can only be reproduced by our models 
if the satellite is accreted along a direct orbit. 



4. Conclusions 

By means of N-body/SPH simulations, we have investigated the 
distribution of orbital eccentricities of stars in thick disks formed 
by heating of a thin stellar disk after minor mergers with 1:10 
mass ratios. The simulated interactions span a range in both 
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Fig. 2. Eccentricity distribution of thick disk stellar particles 
for a 2x(l:10) consecutive minor merger with a gas fraction 
fgas-0-2. As in Fig. 1, stars have been selected in a cylindri- 
cal shell at disk scale lengths, Rj, from 2 < R/Rd < 3, and at 
heights, zo, from 1 < z/zo < 3. The contribution of stars accreted 
from the satellite galaxies is shown in red. 



fgas = 0.2, direct 
Dierickx et al. 2010 
Wilson et al. 2010 
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Fig. 3. Comparison of the eccentricity distributions of stellar 
particles of one of our minor merger mod els with the d i stribu- 
tions found by iDierickx et all (120101) and IWilson et all d201dh 
for stars in the solar neighborhoo d. For comparison, the distribu- 
tion found bv lSales et al.1 (120091) for their "heating" (i.e., minor 
merger) scenario is also shown. In particular, one should note 
the secondary peak at high values of e, which is not seen in our 
simulations of 1 : 10 minor mergers. 



gas fractions of the primary disk and initial orbital parameters. 
The eccentricity distributions have an approximately triangular 
shape, with a peak around e = 0.2-0.35, and a relatively smooth 
decline towards higher eccentricities. The exact distribution de- 
pends on the gas-to-stellar disk mass fraction in the primary 
disk, in the sense that higher gas fractions lead to lower peak 
eccentricity values. It also depends on the orbital configuration 
- direct orbits tend to produce broader eccentricity distributions 
than retrograde encounters. Stars that originated from the satel- 
lite galaxy tend to have relatively high eccentricities after the 
merger is completed (on average e-0.45 to 0.75), similar to those 
recently found in simulations by Sales et al. (2009) but do not 
reach the highest values that these authors found nor show their 
secondary peak at e ~ 0.8. This is likely a consequence of the 
lower satellite mass of the satellites in our simulations, of the dif- 




Fig. 4. Height-dependence of the stellar eccentricity distribution 
in a cylindrical shell with 2 < R/Rd < 3, for a direct (top row) 
and a retrograde (bottom row) 1:10 merger, with f gas = 20% 
. In each panel, the histograms correspond to different heights z 
above and below the galaxy mid-plane as indicated in the legend. 



ferent dark-matter halo profiles used to model the galaxies, and 
our greater range of initial orbital configurations. The absence 
of this secondary peak at high eccentricities, in all our models, 
results in a distributio n of eccentricities that matches qualita - 
tively the observations (IDierickx et alj2010llWilson et alj2010l) . 
Moreover, the increase in the orbital eccentricities of stars in the 
solar neighborhood with vertical distance from the galaxy mid- 
plane found by Diericxk et al. (2010) can also be qualitatively 
reproduced by our models, but only if the merging satellite is ac- 
creted on a direct orbit. No trend is found for retrograde orbits, 
except for stars high above the mid-plane at distances greater 
than 5zo- 

On the basis of these models, we conclude that the observed 
distribution of orbital eccentricities of stars in the Galactic thick 
disk is in qualitative agreement with a scenario where the thick 
disk has been (predominantly?) formed by heating the thin disk 
during the impact of a (or perhaps multiple) minor m erger(s) 
(whose effects are likely cumulative lOu et al ] |20T0ah . As we 
have discussed, the formation of a thick disk in situ, i.e., after 
the merger of gas-rich progenitors in the early universe, cannot 
be ruled out based on our models or other evidence, because its 
observed eccentricity distribution is also qualitatively similar to 
that observed. 

For a satellite accreted o n a direct orbit, our mo dels can also 
reproduce the trend found by IDierickx et alj d2010l) that satellite 
orbits become more eccentric as one moves away from the Milky 
Way plane. We speculate that if minor mergers were the domi- 
nant means of formating the Milky Way thick disk, the major 
contributor should have been merger(s) with satellite(s) on direct 
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orbits. If the thick disk of the Milky Way was formed by external 
heating, other processes may have contributed to its subsequent 
evolution and helped to shape its current characteristics. Radial 
migration dSellwood & Binnevll2002r.lRoskar et all2008albh . for 
example, is particularly efficient if asy mmetries such as bars 
and spiral arms develop in the thin disk ( Debattista et alj|2006t 
iFovle et al.1 120081: iMinchev etaD l2010allbl) and may well have 
played an important role in influencing the kine matic and chem- 
ical p roperties of the early formed thick disk ( ILoebman et alj 
l2010l) . Crucial tests of the viability of processes such as mi- 
nor mergers will be made through detailed comparisons that are 
not restricted to eccentricity data but will involve the kinematic 
properties and metal abundances and ratios of the stars of the 
thick disk. 

The presently available computational facilities ensure that 
the investigation of a large number of parameters and accretion 
histories is feasible. By comparing the signatures produced by 
different accretion events in the kinematics and metallicity of the 
thin/thick disk and halo of the Galaxy, the future observations 
by GAIA and LAMOST, for example, will have a key role in 
understanding the mass accretion history of our Galaxy and its 
impact on the formation and evolution of the thick disk. 
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